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We present a measurement of the helicity of the W boson produced in top quark decays using 
ti decays in the £+jets and dilepton final states selected from a sample of 5.4 fb _1 of collisions 
recorded using the DO detector at the Fermilab Tevatron pp collider. We measure the fractions 
of longitudinal and right-handed W bosons to be fo = 0.669 ± 0.078 (stat.) ± 0.065 (syst.) and 
/+ = 0.023 ± 0.041 (stat.) ± 0.034 (syst.), respectively. This result is consistent at the 98% level 
with the standard model. A measurement with fo fixed to the value from the standard model yields 
/+ = 0.010 ± 0.022 (stat.) ± 0.030 (syst.). 
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I. INTRODUCTION 

o , 

The top quark, which is the heaviest known fundamental particle, was discovered in 1995 0,0 at the Tevatron 
proton-antiproton collider at Fermilab. The dominant top quark production mode at the Tevatron is pp — ¥ tiX. 
Since the time of discovery, over 100 times more data has been collected, providing a large number of ti events 
with which to study the properties of the top quark. In the standard model (SM), the branching ratio for the 
CO | top quark to decay to a W boson and a b quark is > 99.8%. The on-shell W boson from the top quark decay has 

• three possible helicity states depending on the relative direction of the spin and momentum of the W boson. We 
| define the fraction of W bosons produced in these states as fo (longitudinal) , /_ (left-handed) , and /+ (right - 

^ . handed). In the SM, the top quark decays via the V — A charged weak current interaction, which strongly 
suppresses right-handed W bosons. The SM expected values are /o=0.698, /_=0.301, and /+ = 4.1 x 10~ 4 . 
^ . The uncertainties on the SM expectations are ~ (1 — 2)% for fo and /_, and O(10~ 3 ) for / + Q. 

• Here we present a measurement of the W boson helicity fractions fo and /+ and constrain the fraction /_ 
through the unitarity requirement of /_ + / + + fo = 1. Any significant deviation from the SM expectation 

l— ~~ '■ would be an indication of new physics, arising from either a deviation from the expected V — A coupling of the 
I ] tWb vertex or the presence of non-SM events in the data sample. 

The extraction of the W boson helicities is based on the measurement of the angle 9* between the directions 
■^J- ' of the top quark and the down-type fermion (charged lepton or d, s quark) decay product of the W boson in 
lO i the W boson rest frame. The dependence of the distribution of cos#*on the W boson helicity fractions is given 
G\ ■ by 

(N 

w(c) ex 2(1 - c 2 )/ + (1 - c) 2 /- + (1 + cff + (1) 

with c = cos#*. After selection of a ttbar enriched sample the four- momenta of the ti decay products in each 
event are reconstructed as described below, permitting the calculation of cos#*. Once the cosO* distribution 
is measured, the values of fo and /+ are extracted with a binned Poisson likelihood fit to the data. The 
measurement presented here is based on pp collisions at a center-of-mass energy y/s = 1.96 TeV corresponding 
to an integrated luminosity of 5.4 fb^ 1 . 
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II. DATA AND SIMULATION SAMPLES 



This analysis is performed using events collected between April 2002 and June 2009, corresponding to a total 
integrated luminosity of 5.4 fb _1 . Analysis of the Run Ha sample, which totals about 1 fb _1 , was presented in 
Ref. [|| . Here we describe the analysis of the Run lib data sample and then combine our result with the result 
from Ref. [3| when reporting our measurement from the full data sample. 

The Monte Carlo (MC) simulated samples used for modeling the data are generated with ALPGEN [j| interfaced 
to pythia t 6j for parton shower simulation, passed through a detailed detector simulation based on GEANT [7] , 
and reconstructed using the same algorithms as are used for collider data. For the signal (ti) sample, this 
analysis requires MC samples with arbitrary non-standard values for the W helicity fractions, while ALPGEN 
can only produce linear combinations of V — A and V + A tWh couplings. Hence, for this analysis, we use 
samples that are either purely V — A or purely V + A, and use a reweighting procedure to form models of 
arbitrary helicity states. ALPGEN is also used for generating all y+jets processes where V represents the vector 
bosons, pythia is used for generating diboson (WW, WZ, and ZZ) backgrounds in the dilepton channels. 
Background from multijet production is modeled using data. 
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III. ANALYSIS 
A. Event Selection 

For this analysis, the selection is done in two steps. In the first step, a loose initial selection using data 
quality, trigger, object identification, and kinematic criteria is applied to define a sample with the characteristics 
of ti events. Subsequently, a multivariate likelihood discriminant is defined to separate the ti signal from the 
background in the data. We use events in the £+jets and dilepton ti decay channels, which are defined below. 

In the ^+jets decay ti — > W + W~abbb — > iv qq abbb, events contain one charged lepton (where lepton 
here refers to an electron or a muon), at least four jets with two of them being b jets, and significant missing 
transverse energy $t- For the dilepton decay channel, ti — > W + W~ abbb — > Ivl'v'abbb, the signature is two 
leptons of opposite charge, two b jets, and significant $t- 

The main sources of background after the initial selection in the ^+jets channel are VF+jets and multijet 
production; in the dilepton channels they are Z boson and diboson production as well as multijet and VF+jets 
production. The multijet contribution to the ^+jets final states in the initially-selected sample is estimated 
from data following the method described in Ref. In the dilepton channels we model the background due 
to jets being misidentified as isolated leptons using data events where both leptons have the same charge. This 
background originates from multijets events with two jets misidentified as leptons and from IT+jets events with 
one jets misidentified as a lepton. 

TABLE I: The set of variables chosen for use in Lt for the e+jets and /x+ jets channels. The numbers of background 
and ti events in the initially-selected data, as determined from a fit to the L t distribution, are also presented. 





e+jets 


/x+jets 


Events passing initial selection 


1442 


1250 


Variables in best L t 


C 


C 




Ut 


Ut 




K' 

"Tmin 






NN 6avg 


NN ftavg 




xl 


h 










Aplanarity 




N (ti) 


592.6 ± 31.8 


612.7 ± 31.0 


N (IT+jets) 


690.2 ± 21.8 


579.8 ± 18.6 


N (multijet) 


180.3 ± 9.9 


6.5 ± 4.9 



TABLE II: The set of variables chosen for use in Lt for the dilepton channels. The number of background and tt events 
in the initially-selected data, as determined from a fit to the Lt distribution, are also presented. 





e/x 


ee 


nn 


Events passing initial selection 


323 


3275 


5740 


Variables in best L t 


»4,<S,/i,/7ijj m in 
#Tmin#r,NNf,l,TO« 


„4.,<S,7Tljj m i n 

# T ,NNf,i,m« 


xLNNw 


N (ti) 


178.7 ± 15.6 


74.9 ± 10.7 


86.0 ± 13.8 


N (background) 


144.3 ± 14.5 


3200 ± 75 


5654 ± 76 



To separate the ti signal from these sources of background, we define a multivariate likelihood and retain 
only events above a certain threshold in the value of that likelihood. The set of variables used in the likelihood 
and the threshold value are optimized separately for each ti decay channel. The first step in the optimization 
procedure is to identify a set of candidate variables that may be used in the likelihood. We start with twelve 
variabels namely, Aplanarity, Sphericity, Ht, Centrality, K^ min , TO Jjm i n , h, Ac/)(lepton,^T), b jet content 
of the event, |?t or Xz an( i Di- lepton mass ma- We consider all combinations of the above variables to select 
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FIG. 1: (Color online) Best Lt variable for the (a) /x+jets and (b) e+jets channels. The MC is normalized using the 
signal and background fractions returned by the Poisson maximum likelihood fit. The arrows mark the required Lt 
values for events in each channel. 



the optimal set to use for each tt decay channel. For a given combination of variables, the likelihood ratio L t is 
defined as 



«p{E^TNf)f*]} 

L t = r 1 v- 1 — (2) 

exp{E^TK|)f]} + l 

where N vax is the number of input variables used in the likelihood, and (-§)f is the ratio of the parameterized 
signal and background probability density functions. We consider all possible subsets of the above variables 
to be used in L t and scan across all potential selection criteria on L t . For each L t definition and prospective 
selection criterion, we compute the following figure of merit (FOM): 

FOM = h = (3) 

^N s + N B +a% 

where N$ and Nb are the numbers of signal and background events expected to satisfy the L t selection. The 
term erg reflects the uncertainty in the background selection efficiency arising from any mis- modeling of the 
input variables in our MC. 

The signal and background yields in the initially-selected sample for the £+jets channels are listed in Table U 
and for the dilepton channels in Table HU Tables Hill and HVl show the optimal L t cut value for each channel and 
the final number of events in data and the expected numbers of signal and background events after applying 
the L t requirement. Figures [T] and [5] show the distribution of the best likelihood discriminant for each channel, 
where the signal and background contributions are normalized according to the values returned by the fit. 



TABLE III: Expected background and tt yields, and the number of events observed, after the selection on Lt in the 
£+jets decay channels. 





e+jets 


/i+jets 


Best Lt cut 


0.58 


0.29 


Expected tt 


484.4 ± 41.4 


567.2 ± 47.3 


Expected W+jets 


111.7 ± 12.6 


227.7 ± 19.2 


Expected multijet 


58.1 ± 3.9 


4.0 ± 3.1 


Expected total 


656.2 ± 43.4 


798.9 ± 51.2 


Observed 


628 


803 
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FIG. 2: (Color online) Best Lt variable for the (a) efi, (b) ee and (c) fj,fi decay channels. The MC is normalized using 
the signal and background fractions returned by the Poisson maximum likelihood fit to the L t distribution. The arrows 
mark the required L t values for events in each channel. 



TABLE IV: Expected background and tt yields, and the number of events observed, after the selection on Lt in the 
dilcpton decay channels. 







ee 




Best Lt cut 


> 0.28 


> 0.934 


> 0.972 


Expected tt 


186.6 ± 0.4 


44.5 ± 0.3 


43.6 ± 0.3 


Expected Z/-y* ->■ £ + £~ 


N/A 


7.4 ± 1.0 


19.1 ± 1.3 


Expected Z/j* — » tt 


11.2 ± 3.7 


0.8 ± 0.3 


0.35 ± 0.05 


Expected WW 


5.6 ± 1.4 


0.3 ± 0.1 


0.13 ± 0.05 


Expected WZ 


1.5 ± 0.5 


0.28 ± 0.04 


0.16 ± 0.01 


Expected ZZ 


1.0 ± 0.5 


0.34 ± 0.04 


0.57 ± 0.04 


Expected misidentified jets 


15.9 ± 3.1 


0.54 ± 0.48 


3.7 ± 2.5 


Expected total 


221.7 ± 5.1 


54.2 ± 1.2 


67.7 ± 3.9 


Observed 


193 


58 


68 



B. Templates 

After the final event selection if performed, cos#* is calculated for each event by using the reconstructed top 
quark and W boson four-momenta. In the i+jets decay channel, the four-momenta are reconstructed using a 
kincmctic fit with the constraints: (i) two jets should give the invariant mass of the W boson (80.4 GeV/c 2 ), (ii) 
one lepton and the $t should give the invariant mass of the W boson, and (iii) the mass of the reconstructed 
top and anti-top quark should be 172.5 GeV/c 2 . In the £+jets decay channel, the hadronic W boson decay 
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FIG. 3: (Color online) Comparison of the cos#* distribution in Run lib data and the global best-fit model (solid line) 
and the SM (dashed line) for (a) leptonic W boson decays in ^+jets events, (b) hadronic W boson decays in £+jets 
events, and (c) dilepton events. 



from the top quark in the event also contains information about the helicity of that W boson. Since we do not 
distinguish between jets formed from up-type and down-type quarks, we choose one of the W boson daughter 
jets at random as the basis for the calculation. With this choice, left-handed and right-handed W bosons have 
identical |cos#*| distributions, but we can distinguish either of those states from longitudinal W bosons, thereby 
improving the precision of our measurement. 



C. Model-independent W Helicity Fit 



The W boson helicity fractions are extracted by computing a binned Poisson likelihood L(fo, /+) with the 
distribution of cos8* in the data to be consistent with the sum of signal and background templates. The 
likelihood is a function of the W boson helicity fractions /o and / + . While performing the fit, both /o and /+ 
are allowed to float freely, and the measured W helicity fractions correspond to those leading to the highest 
likelihood value. The comparison between the best-fit model and the data is shown in Fig. [3l 
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D. Systematics 

Systematic uncertainties are evaluated using simulated event ensembles in which both changes in the back- 
ground yield and changes in the shape of the cosO* templates in signal and background are considered. The 
simulated samples from which the events are drawn can be either our nominal samples or samples in which 
the systematic effect under study has been shifted away from the nominal value. In general, the systematic 
uncertainties assigned to fo and /+ are determined by taking an average of the absolute values of the differences 
in the average fit output values between the nominal and shifted V — A and V + A samples. 

TABLE V: Summary of the systematic uncertainties on /+ and fo- 



Source 


Uncertainty (/+) 


Uncertainty (fo) 


Jet energy scale 


0.007 


0.009 


Jet energy resolution 


0.004 


0.009 


Jet ID 


0.004 


0.004 


Top quark mass 


0.011 


0.009 


Template statistics 


0.012 


0.023 


ti model 


0.022 


0.033 


Background model 


0.006 


0.017 


Heavy flavor fraction 


0.011 


0.026 


b fragmentation 


0.000 


0.001 


PDF 


0.000 


0.000 


Analysis consistency 


0.004 


0.006 


Muon ID 


0.003 


0.021 


Muon trigger 


0.004 


0.020 


Total 


0.032 


0.060 



E. Results 

Applying the model independent fit to the Run lib data, we find 

fo = 0.739 ± 0.091 (stat.)± 0.060 (syst.) (4) 
/+ = -0.002 ± 0.045 (stat.)± 0.032 (syst.). (5) 

The 68% and 95% C.L. contours in the (f+,fo) plane are shown in Fig. @] (a). Finally, we perform fits in 
which one of the two helicity fractions is fixed to its SM value. Constraining fo, we find 

/+ = 0.014 ± 0.025 ± (stat.) ± 0.028(syst.), (6) 

We also constrain / + and measure fo, finding 

fo = 0.735 ± 0.051 (stat.) ± 0.051(syst.). (7) 

F. Combination with Our Previous Measurement 

To combine this result with the previous measurement from Ref. we repeat the maximum likelihood fit 
with the earlier and current data samples and their respective MC models, treating them as separate channels 
in the fit. This is equivalent to multiplying the two-dimensional likelihood distributions in fo and /+ corre- 
sponding to the two data sets. We determine the systematic uncertainty on the combined result by treating 
most uncertainties as correlated (the exception is template statistics) and propagating the uncertainties to the 
combined result. The results are presented in Table fVTl 

The combined result for the entire 5.4 fb _1 sample is 



fo = 0.669 ± 0.078 (stat.) ± 0.065 (syst.) 
/+ = 0.023 ±0.041 (stat.) ±0.034 (syst.). 



(8) 
(9) 
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TABLE VI: Summary of the combined systematic uncertainties on /+ and /o for Run Ha and Run lib. 



Source 


Uncertainty (/+) 


Uncertainty (/o) 


Jet energy scale 


0.009 


0.010 


Jet energy resolution 


0.004 


0.008 


Jet ID 


0.005 


0.007 


Top mass 


0.012 


0.009 


Template statistics 


0.011 


0.021 


tt model 


0.024 


0.039 


Background model 


0.008 


0.023 


Heavy flavor fraction 


0.010 


0.022 


b fragmentation 


0.002 


0.004 


PDF 


0.000 


0.001 


Analysis consistency 


0.004 


0.006 


Muon ID 


0.002 


0.017 


Muon trigger 


0.003 


0.024 


Total 


0.034 


0.065 



(b) 

D0, L = 5.4 fb" 1 

• Best-fit value 

* SM value 



-0.2 0.2 0.4 0.6 0.8 1 1.2 ""-0.2 0.2 0.4 0.6 0.8 1 1.2 



FIG. 4: Result of the model-independent W boson helicity fit for (a) the Run lib data sample and (b) the combined 
Run Ha and Run lib data sample. In both plots, the ellipses indicate the 68% and 95% C.L. contours, the dot shows 
the best-fit value, the triangle corresponds to the physically allowed region where fo + f+ < 1, and the star marks the 
expectation from the SM. 

The 68% and 95% C.L. contours in the (/+, /o) plane are shown in Fig. [4] (b). The probability of observing a 
greater deviation from the SM due to random chance is 83% when only statistical uncertainties are considered 
and 98% when systematic uncertainties are included. 
Constraining /o to be 0.7, we find 

/+ = 0.010 ± 0.022 (stat.) ± 0.030 (syst.) (10) 




Constraining /+ to gives 

fo = 0.708 ± 0.044 (stat.) ± 0.048 (syst.) 



(11) 
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IV. CONCLUSION 



We have measured the helicity of W bosons arising from top quark decay in tt events using both the £+jets 
and dilepton decay channels and find 

f = 0.669 ± 0.102[±0.078 (stat.) ± 0.065 (syst.)], 

(12) 

/+ = 0.023 ± 0.053[±0.041 (stat.) ± 0.034 (syst.)]. 

in a model-independent fit. The consistency of this measurement with the SM values /o = 0.698, /+ = 3.6 x 10 -4 
is 98%. Therefore, we report no evidence for new physics at the tWb decay vertex. 
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